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The properties of the thermal Abelian monopoles are studied in the deconfinement phase of 
the SU(2) gluodynamics. To remove effects of Gribov copies the simulated annealing algorithm 
is applied to fix the maximally Abelian gauge. To study monopole profile we complete the first 
computations of excess of the nonabelian action density as a function of the distance from the 
center of the thermal Abelian monopole. We have found that starting from the distances w 2 lattice 
spacings the chromoelectric and chromomagnetic action densities created by monopole are equal to 
each other, from what we draw a conclusion that monopole is a dyon. Furthermore, we find that the 
chromoelectric and chromomagnetic fields decrease exponentially with increasing distance. These 
findings were confirmed for different temperatures in the range T/T c € (1.5,4.8). 
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One of the hypotheses which have been put for- 
ward in the recent past is that the quark-gluon plasma 
(QGP) properties might be dominated by a magnetic 
component [IHl]- The monopoles or center vortices 
might be responsible for unexpected properties of the 
hadron matter at T > T c : on one hand it is well known 
from lattice results that the equation of state is close 
to that of an ideal gas, on the other hand the very 
low viscosity to entropy ratio tells that it is an ideal 
liquid. 

In Ref. such magnetic component has been re- 
lated to thermal Abelian monopoles evaporating from 
the magnetic condensate which is believed to induce 
color confinement at low temperatures. Moreover it 
has been proposed to detect such thermal monopoles 
in finite temperature lattice QCD simulations, by 
identifying them with monopole currents having a 
non-trivial wrapping in the Euclidean temporal direc- 
tion US!. 

The way one can study the monopoles properties 
on the lattice is via an Abelian projection after fix- 
ing the maximally Abelian gauge (MAG) This 
gauge as well as the properties of the monopole clus- 
ters has been investigated in a numerous papers both 
at zero and nonzero temperature (see for extensive 
list of references, e.g. |9|). The evidence was found 
that the nonperturbative properties of the gluody- 
namics such as confinement, deconfining transition, 
chiral symmetry breaking, etc. are closely related 
to the Abelian monopoles defined in MAG. This was 
called a monopole dominance. 

First numerical investigations of the wrapping 
monopole trajectories were performed long ago in 
Refs. [J] and [5|. A more systematic study of the ther- 
mal monopoles in SU(2) Yang-Mills theory at high 
temperature has been performed in Refs. |R]| - [l2l |. In 



particular, it was found in [Io| that the density of 
monopoles is independent of the lattice spacing, as it 
should be for a physical quantity. 

In paper fill ] very interesting properties of the ther- 
mal Abelian monopoles were found. The authors 
measured the excess of chromoelectric and chromo- 
magnetic action density on the surface of the lattice 
(hyper-)cubes with monopoles inside. The depen- 
dence of the excess on the distance from the monopole 
center was determined through the variation of the 
lattice spacing a (similar investigation at T = was 
made in [Hj]). As a result it was found that with good 
accuracy the chromomagnetic and chromoelectric ac- 
tion densities created by monopole have the following 
behavior: H 2 (r),E 2 (r) = ojj^/r 4 . The coefficients 
an and a e turned out to be equal to each other with a 
very good accuracy, from what the authors concluded 
that monopole is a dyon. It is worth to note that there 
were other works in the past where dyonic properties 
of the monopoles were observed fl4Kla | . 

The drawback of the study undertaken in |Il[ is that 
all results were obtained at the ultraviolet cutoff scale 
and were thus subjected to both lattice discretization 
errors and ultraviolet divergences. In view of the im- 
portance of the findings of [Tl| in this paper we are go- 
ing to study the chromoelectric and chromomagnetic 
fields created by monopole and to check whether the 
observed behavior is correct or is just a lattice arti- 
fact, i.e. artifact of the ultraviolet cut off. To accom- 
plish this check we will measure the chromoelectric 
and chromomagnetic fields at various distances from 
the monopole center. 

In this paper we study the SU(2) lattice gauge the- 
ory with the standard Wilson action 
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where p = & n( i 50 is a bare coupling constant. 
The link variables U x ^ £ SU{2) transform under 
gauge transformations g x as follows: 



U 



Xfi 



g x £ SU(2) . (1) 



Our calculations were performed on the asymmetric 
lattices with lattice volume V = LtL s , where Lt, s is 
the number of sites in the time (space) direction. The 
temperature T is given by 
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where a is the lattice spacing. 

The MAG is fixed by finding an extremum of the 
gauge functional 



Fu{g)= —Y- 
u ™> AV ^ 2 



Tr 



(3) 



with respect to gauge transformations g x . We apply 
the simulated annealing (SA) algorithm which proved 
to be very efficient for this gauge I17ll as well as for 
other gau ges such as center gauges [18( and Landau 
gauge [13. To further decrease the Gribov copy ef- 
fects we generated 10 Gribov copies starting every 
time gauge fixing procedure from a randomly selected 
gauge copy of the original Monte Carlo configuration. 

In Table U we provide the information about the 
gauge field ensembles used in our study. 





a[fm] 


L t 


L a 


T/T c 


^meas 


2.43 


0.108 


4 


32 


1.5 


1000 


2.5115 


0.081 


4 


28 


2.0 


400 


2.635 


0.054 


4 


36 


3.0 


500 


2.80 


0.034 


4 


48 


4.8 


400 



TABLE I: Values of /3, lattice sizes, temperatures, num- 
ber of measurements and number of gauge copies used 
throughout this paper. To fix the scale we take yfa = 440 
MeV. 



The chromomagnetic action density at a site x is 
defined as 



(4) 



The sum is taken over all spatial plaquettes P s which 
contain the lattice site x. In the continuum limit this 
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expression is proportional to ~ T^G^ 
Tr(Hf + H\ + -fff ). So, this expression can be taken 
as a measure of the chromomagnetic action Tr(H 2 ) at 
the site x. 

Analogously, for the chromoelectric action density 
at a site x we take 



Se{x) 
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Here the sum is taken over all time-like plaquettes 
Pt which contain the site x. In the continuum limit 
this expression is proportional to ~ Tt(Gqj + G§ 2 1 
G n3 ) 



os j — Tr(E{ + E£ + £3) and thus it can be taken as 
a measure of the chromoelectric action Tr(E 2 ). Note 
that our definitions for Sm,e(%) differ from those used 
in Ref. [ll|. Although definitions of [ll| were natural 
for the surface of a cube with monopole our definitions 
are more suitable for measurements at some distance 
from such cube. 

Since we are studying the fields created by a 
monopole we should subtract the vacuum fluctua- 
tions of the chromomagnetic and chromoelectric ac- 
tions from the equations ([3]), We define the excess 
of the action density as 



(SSm,e(cI)} = {Sm,e(x)) - (S m .e) 



(0) 



where (...) means ensemble average, d is the distance 
from a monopole and bar means averaging over all 
wrapped monopoles and all sites x at the distance d 
from monopole centers. 

The monopole currents and their wrapping numbers 
are defined in a standard way (see e.g. [10|h Moving 
along wrapped monopole clusters on a dual lattice we 
detect all 3-dimensional cubes in all time slices on the 
original lattice which contain monopoles correspond- 
ing to J4 component of the magnetic current. Having 
detected all such 3-dimensional cubes we do the mea- 
surements of the chromomagnetic and chromoelectric 
action densities (6SM,E(d)) at various distances from 
a given monopole and then we take the average over 
all thermal monopoles found on the lattice. 

Now let us consider a three dimensional cube with 
the monopole belonging to a wrapped cluster. Be- 
low it will be assumed that the monopole is located 
in the center of this cube and we take the center as 
a coordinate origin. Actually, one cannot assert that 
the monopole is exactly located at the center of the 
cube. However, since we take an average over all 
monopoles, this approximation can be considered as 
a good one. We have measured the action densities 
(5Sm e{A)) at various distances d from the centers of 
the cubes with monopole. The distance was defined as 
a length of the vector d = {ni±l/2, ri2±l/2, ri3±l/2} 
from the coordinate origin to lattice site x under 
consideration. We present results of measurements 
for d — ^{m,m,m},rn =1,3 and d = ^{1,1,3}, 
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FIG. 1: (SSm) defined in eq.© as function of the di- 
mensionless distance rT from the monopole center at the 
temperatures T/T c = 1.5, 2.0, 3.0, 4.8. 

±{1,3,3}, ±{1,1,5}, |{1,3,5}. For T/T c = 1.5 ad- 
ditionally results for d = |{1, 5, 5}, 1, 7} are pre- 
sented. For longer distances the statistical errors were 
too large. 

In Figures Q] and [2] results are shown for the (SSm) 
and (SSe) as functions of the dimensionless distance 
rT = d/4. From these Figures one clearly sees that 
at least at large distances both (SSm) and (SSe) de- 
crease with distance in agreement with exponential 
fall-off ~ exp (— 2M m<e r). The dependence 1/r 4 found 
in [ll| is ruled out. We do not have enough data points 
to determine the pre-exponential function by fitting. 
Respectively, it is rather difficult to find the parame- 
ters M mte with a good accuracy. In this paper we just 
make rough estimation of these parameters fitting the 
last 3 data points for (SSm) to the exponential fall-off 
with constant prefactor. We get the following results: 
M m /T = 3.5(2), 4.0(4), 4.3(2), 3.7(8) for the tempera- 
tures T/T c = 1.5,2.0,3.0,4.8, respectively. 

Looking at Figures [T] [2] one can see that the data 
lie on the smooth curves. This means that the data 
obey rotational invariance, since the data at different 
distances were measured in different directions. More- 
over, vectors d = ^{3,3,3} and d = ±{1,1,5} have 
equal length and one can check the rotational invari- 
ance directly. Indeed, we find for these two vectors 
consistent results with deviations within 2a interval. 
In all Figures we show averaged data for these vectors 
d. 

At large enough distances ( beginning from the 
distance d = 2.18 ) the monopole chromomagnetic 
and chromoelectric action density seem to be equal 
to each other. To demonstrate this important prop- 
erty we plot the ratio (SS M )/(SS E ) = H 2 /E 2 in Fig- 
ure |3l From this Figure we see that within the error 
bars at distances rT > 0.5 the ratio (SSm) / (SSe) is 
compatible with 1. From this observation one can 



FIG. 2: Same as in Figure [T] for (SSe). 
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FIG. 3: The ratio (5S M )/(SS E ) = H 2 /E 2 as function of 
the distance rT from the monopole center at temperatures 
T/T c = 1.5,2.0,3.0,4.8 The data sets at T/T c = 2.0 and 
3.0 are shifted along horizontal axis to improve readability 
of the figure. 

draw a conclusion: at least at large enough distances 
H 2 (r) — E 2 (r). This implies that monopoles carry 
both chromoelectric and chromomagnetic charges and 
they are equal. So monopoles are selfdual dyons. 
These statements are the main results of this paper. 

To understand the reason of lack of selfduality at 
small distances let us look at Figure |3] For all tem- 
peratures we observe similar behavior. At the distance 
d = 0.87 in lattice units the ratio (SSm) /(SSe) ~ 2 
for all temperatures. At the distance d = 1.66 the ra- 
tio (SSm) /(SSe) ~ 1-3. At larger distances the ratio 
is compatible with unity. We believe that deviation of 
the ratio from one at small distances can be explained 
by the discretization effects. Notice that our defini- 
tion of the chromomagnetic and chromoelectric action 
densities at site x is nonlocal involving all plaquettes 
in respective planes which own site x. This nonlocal- 
ity is different for two action densities: it is purely 
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spatial for the chromomagnetic action density and is 
both spatial and temporal for the chromoelectric one. 
Thus at distances of order of one lattice spacing we, 
evidently, measure the fields taken at different points 
and different distances. These our arguments should 
be checked by computations with smaller lattice spac- 
ing, i.e. with Lt > 4. This will be done in a forth- 
coming paper. In that paper we will also present our 
data on the density and interactions of the thermal 
Abelian monopoles. 

It is clear that results and conclusions of Ref. [Tlj 
where (SSm) and {6Se} were measured at the near- 
est possible distance to the monopole center are sub- 
jected to same discretization effects as discussed above 
for our data at small distances. While the distance 
dependence I/r 4 found in Ref. [ll[ is an ultraviolet 
divergence effect. 

Thus we established that the Abelian thermal 
monopoles carry both chromoelectric and chromo- 



magnetic charges and they are equal. So monopoles 
are selfdual dyons. Furthermore respective action 
densities are screened. We believe that these re- 
sults are important for understanding QCD in the 
quark-gluon plasma phase since many recent theoret- 
ical models of this phase include monopoles as an im- 
portant ingredient. These are model of Ref. 0, Q 
based on competition between magnetic and electric 
quasiparticles, dyon [20| and caloron [2l[ models. 
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